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Abstract 

A diagrammatic approach based on thermal Green function to su- 
perconducting fluctuation transport is reviewed keeping consistency 
with Ginzburg-Landau theory. The correct expression of the heat 
current vertex for Cooper pairs is clarified via Jonson-Mahan formula 
and Ward identities. 

1 Introduction 

In this Note, I, a diagrammer, try to supply you, beginners, with some de- 
tailed calculations usually skipped in the references written by professionals^ 
I hope that you can save time and efforts by this note and reserve enough 
energy to go farther. 

This note is intended to be a primei^l to the textbook: Theory of Fluc- 
tuations in Superconductors [2|, [3] . It is the first of the series consisting of 
three notes. This first note mainly discusses the transport in the absence of 
magnetic field. The second and the third mainly discuss the transport under 
magnetic field. 



*Departnient of Physics, Kyushu University, Fukuoka 812-8581, Japan 
^ I shall remark on some mistakes made by professionals in several footnotes. Such 
notes might be helpful for you to avoid unnecessary struggle with the mistakes. 

^ In the discussion of superconducting fluctuation transport (T > Tc) I only consider the 
Aslamazov-Larkin (AL) contribution and leave the consideration of the Maki- Thompson 
(MX) and the density-of-states (DOS) contributions to the textbook . The AL trans- 
port is equivalent to the time-dependent Ginzburg-Landau (GL) transport. 
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I expect that you have some knowledge of the Feynman-diagram technique 
at finite temperature and the hnear response calculation via the thermal 
Green function. 1 will skip detailed explanations of these in the following. 

We use the unit of h = ks = c = 1. 

2 Hamiltonian 

In grand canonical formalism we use K 

K^H-i^N, (1) 

instead of the Hamiltonian H where N and /i are the number and chemical 
potential of electrons. We decompose K as 

K = Ko + V, (2) 

where 

p 

is the kinetic energy measured from the chemical potential jj, with 
and 

q 

is the local attractive interaction with g being a positive constant. Here 
represents the creation operator of an up-spin electron with momentum p, 
bp represents the annihilation operator of a down-spin electron and so on. 
The creation and annihilation operators of Cooper Pairs are 

p p 

The (imaginary) time dependence of the operator A is given as 

A{t) = e^Me-^^ (7) 

or 

§^A{r)^[K,A{r)]. (8) 
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3 Electron Propagator 

The propagator for electrons G{p, r) is defined as 

G(p,r) = -(T4ap(r)at,}>, (9) 

where (X^ is the expectation value of X using the density matrix defined 
by K. In the following, throughout the series of three notes, the Zeeman 
splitting is neglected so that 

-{T^{b^{r)bl}) = G{p,r). (10) 

The operation of the (imaginary) time ordering is explicitly expressed 



as 



G{p,t) = -e{T){a^{T)al) +e{-T){ala^{r)), (11) 

for electrons where 6{x) = 1 for x > and 6{x) = for x < 0. Since the 
time evolution of the free electron is given a^ 

ap(r) = e-«p^ap, (12) 

the propagator for free electrons Go(p, t) become^ 

Go(p, r) = e-«p^{ - e{r) [l - /(^p)] + e{-r)fi^^)}, (13) 

where f{x) is the Fermi distribution function 

/(^) = (14) 

The full propagator G(p, r) is expressed as the superposition of this-type of 
function: 

/oo 
dep{e)e-^^{-eiT)[l-fie)]+ei-T)fie)}. (15) 
'OO 

Here p(e) is the spectral function. 

The Fourier transform of the propagator is defined as 

^ CO 

^(P'^) = ^ E G(p,^e„)e--'^^ (16) 



3 alir) = e«p-a^. 



" Qpuj^ = 1 — a^pQp. For free electrons (ajjap)^ = f{S,p) where (^)q is the expectation 
value of X using the density matrix defined by Kq. 
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G{^,ten)= / drG'(p,r)e^^"^ (17) 
Jo 

where = {2n + l)7rT with n being integer. Here T is the temperature and 
/3 is the inverse temperature /3 = 1/T. 
The free propagator is transformed as 

Go(p, ^e^) = - f dr(ap(r)ai>e^^"- = " [l " /(^p)] rfre^^^--?-)^ (18) 

JQ Jo 

Using e*^"^ = —1, we obtain 

Goip^tsn) = . \ . (19) 

^^n Sp 

The full propagator G{p, iSn) is expressed as the superposition of this-type 
of function: 

G(p,^£„)= r de-^^. (20) 

The retarded and advanced propagators, G^ijp.e) and G^(p,e), are ob- 
tained from the thermal propagator G'(p, 

G^(p,£) =G(p,£ + z5), G^(p,e) =G(p,£-25), (21) 

where 5 = +0. The spectral function p{e) is relatecH to the imaginary part 
of the retarded propagator as 

ImG^(p,e) = -7rp(e). (22) 

The full and free propagators are related by the Dyson equation 

^(p, i£nY^ = Go{p, iEn)'^ - S(p, «£„), (23) 

where S(p, z£:„) is the self-energy. 

4 Cooper-Pair Propagator 

The effective interaction for Cooper pairs in the ladder approximation L(q, iUrr^ 
is given by 



-L(q, ium) = g + gUiq, iuJm)g + 5'n(q, iujm)g'n{q, iuJm)g + 
^ 9 

1 - 5'n(q, iuj^y 



(24) 
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We have employed the identity 



I-KO^X)^ 



X + iS X 

where P denotes that the principal value should be taken when it is integrated. 
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Figure 1: Effective interaction for Cooper pairs: Tlie wavy and broken lines 
represent —L and —{—g). The solid line with arrow represents Gq. All the 
Feynman diagrams in this note are drawn by JaxoDraw. 



wher^ 

n(q, ium) = T^^ G'o(p + ^, 2e„)G'o(-P + ^, i^^m - iSn), (25) 

n p 

and ujrn = 2m7rT with m being integer. 

The propagator for Cooper pairs -D(q, r) is defined as 

D{q,r) = -{Tr{P^{r)Pl}), (26) 

and its Fourier transform in the ladder approximation D{c[,iUm) is given by 



-D{ci, ium) = n(q, ium) + n(q, iUm)gll{ci, iUm) + 
n(q, iUm) 



(27) 



1 - 5'n(q, iuJm)' 
The relatior0 between D{q^,iujm) and L(q^,iu„i) is 

L(q, ium) = -g + g^Diq., i^m)- (28 



^ Since the particle-particle pair propagator n(q, icj„i) depends on the shape of the 
region of p-summation, we have chosen the symmetric configuration in p-space as (j25p . 
On the other hand, the particle-hole pair propagator is free from the choice. See Fukuyama, 
Hasegawa and Narikiyo: J. Phys. Soc. Japan 60, 2013 (1991) for these discussions. The 
imaginary part of P^{q, to) proportional to derived in Sergeev, Reizer and Mitin: Phys. 
Rev. B 66, 104504 (2002) is absent in the symmetric choice. 

^ In the forthcoming discussion employing the Ward identity the inverse of the propa- 
gator -D(q, iuj„i)^^ plays the central role. In the limit of T — > Tc, q ^ and w„i — !> 0, 

D{q,iuj^y^ = L{q,iuj^y^ ■ =. L(q,ia;™)~^ • ^ = L(q,ia;™)"^ • 5^ 

n(q,icj„i) n(0, 0) 

where 1-50(0,0) ^0{T-Tc) and i(q, itj„)-i - 0(r - T^, q^, a;,„) as 

L(q, lujm)-' = - [g-^ - n(o, 0)] - [n(o, 0) - n(q, luj^)] . 

Thus in this limit 

£'A(q, it^m)"^ Ni(q,iw„)"\ 
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In terms of the propagator of the order-parameter D/^[c\,iUm) this relation 
is rewritten as 

L(q, ium) = -g + DA{q, iuJm), (29) 
where the propagator is introduced as 

DA{q,r) = -{T^{A^{r)A\}), (30) 

with 

A{ = g-Pl, A^ = g- P^. (31) 

The retarded and advanced propagators, D^{q,ijj) and D^{q,co), are ob- 
tained from the thermal propagator D{q,iUm) as 

D^{q,uj) = D{q,uj + i5), D^{q,uj) = D{q,UJ - i5). (32) 

In the ladder approximation, the full and free propagators, D{q, iUm) and 
DQ{q,iuJrn), are related by the Dyson equation as 

D{q, iUm)'^ = Do{q, iUm)'^ - S(q, iUm), (33) 

where 

Do{q,iUm) = -'n{q,ioJm), T.{q,iujm) = -g- (34) 

In the limits of low energy and long wavelength the retarded effective inter- 
action L^{q, to) = L(q, to + i5) obtained from i^ 

^) = - .L, ^.2 2 • ' (35) 

^(0) e + CoQ - ^^^0 

where 

T T — T 

e = ln4=.i— ^, (36) 

and A^(0) is the density of states per spin at the chemical potential. This 
effective interaction has Ornstein-Zernike form in space-direction and Debye 
form in time direction. The functional form of ( 135|) holds beyond the ladder 
approximation so that we regard e, and tq as phenomenological parameters 
in the following. 



^ See §6.2 in j3] for the microscopic calculation in the ladder approximation. In §6.4 of 
[3] the renormalization due to impurity scattering is also discussed. 
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5 Boltzmann Transport: Relaxation- Time Ap- 
proximation 

In the relaxation-time approximation of the Boltzmann transport^ the ex- 
pectation value of the charge current and the heat current J*^ are given 
aS 

Jl = (e^'o;) = 2e ^ v^Qp, (37) 
p 

and 

= {i^v.)='^Y.^pV,3^. (38) 
p 

where is the deviation of the distribution functioiJ"! from the equilibrium 
distribution function f{i-p) defined by f|T4|) and e is the charge of an elec- 
tron (e < 0). Throughout the series of three notes the Zeeman splitting 
is neglected so that the spin degrees of freedom is accounted only by the 
degeneracy factor 2. 

If a static electric field E = {E^, 0, 0) is applied, the resulting deviation is 

5p = eE.t;.r(-^^), (39) 
J: = 2e'E. E - ^) - (40) 



so that 



In the same manner 



p ^ ^ 

If a static magnetic field H is applied additionally, the resulting is 



(-^^)fE--(HxE)l, (42) 



+ {uj^tY V dip 



^See, for example, Ziman: Principles of the Theory of Solids, 2nd edition (Cambridge 
Univ. Press, Cambridge, 1972) for the Boltzmann transport of electrons. The Boltzmann 
transport of Cooper pairs is discussed in §3.7 of [S . 

10 Here J- = ( J|, Jf ), JQ = ( jQ, jQ, jQ) and Vp = {v,,Vy, v,). 

11 Since we describe the distribution function of quasi-particles by the Boltzmann equa- 
tion, and v-p represent the energy and the velocity of an quasi-particle with a renormal- 
izedmassm* within the discussion of the Boltzmann transport. Namely, ^p = vf-{\v\~Pf) 
and Vp = p/m* where pf is the Fermi momentum and vp — pp/m* . To keep the notation 
simple, we do not distinguish m* and m in the following. 
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where Uc = |eH|/m. 

If a static temperature gradient VT is applied additionally, its effect is 
taken into account as 



eE' = eE - 



VT 



so tlia 



and 



jQ _ 2e V ^ f 



dm 



E'- — (H X E') 

m 



(43) 



(44) 



PT 

E'-— (HxE') . (45) 



6 Boltzmann Transport: Exact Formula 

The Boltzmann equation in linear order of E for an interacting electron 
system under static electromagnetic field is 

eE.Vp^ + e(vpXH).^ = Cp, (46) 

where the collision term Cp is given as 

= 5Z {^PP'V - ^P'P^p} = (^t7^)pp/V> (47) 



with 



and 



(^t7')pp, = :;^W-Cpp', (48) 



-^E^P'P- (49) 



'P 

P 



The treatment of the magnetic field in Reizer and Sergeev: Phys. Rev. B 61, 7340 
(2000) is not compatible with the Boltzmann transport. That in Sergeev, Reizer and 
Mitin: Phys. Rev. B 77, 064501 (2008) is not compatible with the GL transport. In these 
works the gauge invariance is not satisfied in ordinary manner where the momentum is 
combined with the vector potential and the energy with the scalar potential. While the 
direction of the electron motion is affected, the value of the energy is not affected by the 
magnetic field. Some criticisms against these works have been made by Serbyn, Skvortsov 
and Varlamov: arXiv:101 2.4316 However, Sergeev, Reizer and Mitin: . arXiv :1101.4186 
still claim their validity. 
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Here Cp'p is the transition rate from p to p' and has a symmetry Cp'p = Cpp' . 

This hnear Boltzmann equation is formally solved exactly!^ as follows. 
By introducing the matrix A 



^PP' = (nr')pp, -e(v. 
6]) is written in the form 



H 5p,p', 



P' 



(50) 



(51) 



If we get the inverse matrix A ^, the deviation of the distribution is known 
as 



(52) 



Using the identity for operators P and Q 

(P -Qy^ = P-' + P-^QP-^ + P-'QP-^QP-^ + ■ • •, (53) 
the inverse matrix is obtained as an expansion in terms of the magnetic field 

d 



Vp = (-tr)p,p + e J2 inr) p,p^ (vp, X ^) -H (r,)p^^ 
Pi L 



Pl P2 











[K"9pj 


H 







Employing f l52|) and fl54l) we can calculate 



and the conductivity tensor a^'^ is introduceco as 



Thus the conductivity tensor is expressed as 



p' P 



(^-)p.p + -- 

(54) 
(55) 

(56) 
(57) 



See, for example, Kotliar, Sengupta and Varma: Phys. Rev. B 53, 3573 (1996). 
Here = (J|, J^, Jf), E = {E^,Ey,E^) and Vp = (Wp,w^,Up) with = x,y,z. 
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and this expression is microscopically derived on the basis of the Fermi-liquid 
theory!^ The conductivity in the absence of magnetic field is obtained as 

a« = 2e^5:5:„;,(n,),,„«5(-^). (58) 



7 Charge and Heat Currents 

The charge current of electrons, created by 'ipli'x.) and annihilated by V^(t(x), 
is 

j^W = KW(^^'^W) - (vV^i(x))v^.(x)], (59) 

and satisfies the conservation law 

p(x) + V-r(x)=0, (60) 

where 



is the charge density. Using the four- divergence (IBUl) is written as 

E^^'^(^) = 0> (62) 

11=0 ^ 

where x = (t, x) = {xo^i,X2, x^) with /i = 0, 1, 2, 3. 
The energy currenio of electrons is 

"(^) = -^E [(^^i(^))^'^(^) +^i(^)(^^'^(^))]' (63) 



The microscopic derivation of the conductivity is done by Ehashberg: Sov. Phys. 
JETP 14, 886 (1962). This work is extended to the case of the Hall conductivity by 
Kohno and Yamada: Prog. Theor. Phys. 80, 623 (1988) and to the magneto-conductivity 
by Kontani: Phys. Rev. B 64, 054413 (2001). 

Via the energy-momentum tensor the energy current is given as 

where = 1, 2, 3. The Lagrangian density C can be replaced by its kinetic part 

because its interaction part Cint does not contain the derivative of ipl. and il^a- See, for 
example, Langer: Phys. Rev. 128, 110 (1962) or §10.3 in [5]. 



10 



and satisfies the conservation law 

/i(x) + V ■ u(x) = 0, (64) 

where h{x.) is the Hamihonian density. 

According to the thermodynamic relation dQ = dU—fiN, the heat current 

is 

j«(x) = u(x) - ^j^(x), (65) 

e 

and satisfies the conservation law 

3 

dx 



E9^#(^) = 0' (66) 

fi=0 ^ 



where 

Jo^(x) = Mx) - ^p(x). (67) 

8 Kubo Formula: Linear Response to Elec- 
tric Field 

We observe the expectation value of the charge current J"^ caused by the 
external electric field E as 



j;(k,a;) = 5^a^,(k,u;)E,(k,u;). (6^ 



The conductivity tensor a is given by the Kubo formulaLj (linear response 
theory) 



<^%{k,zuj,)= [ dTe'-^^{Tr{j;{k;r)j:{k = 0)}), (69) 

and 

a,,{k, u) = ^ [$^,(k, u + zS)- $^,(k, z5)] . (70) 

Here 

j:ik;T) = e^^fJk)e-^\ (71) 



^"^ I follow Fukuyama, Ebisawa and Wada: Prog. Theor. Phys. 42, 494 (1969) concern- 
ing the use of the Kubo formula. The sign of is opposite to Q^^ in [3 and [AGD] and 
to Kf,^ in "4" and [FW]. 

[AGD] = Abrikosov, Gorkov and Dzyaloshinski: Methods of Quantum Field Theory in 
Statistical Physics (Dover, New York, 1975). 

[FW] = Fetter and Walecka: Quantum Theory of Many- Particle Systems (McGraw-Hill, 
New York, 1971). 
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and 

j^(k) = / £xe-'^''i'{x). (72) 



The Fourier transform of ( 159|) results in 

k 



p 

= e 5^ ^ (vp + Vp+k) (aj,ap+k + &|,6p+k) ' (73) 



2 

p 



and in the hmit of k — t- 

f(k = 0) = e5^vp(at,ap + 6i6p) =e5^5^vpct,^cp., (74) 



P cr 



where c|,g. represents the creation operator of an electron with momentum p 
and spin a and Cp„ represents the annihilation operator (a =t, I) introduced 

by 

(x) = e-^^'-^ct ^, M^) = E ^^^"'^P- C^S) 

p p 

In the presence of the temperature gradient both E and VT are the causes 
of the observed charge current and heat current 3^ as [S] 



3^ \ f a a \ f E 



(76) 



where k is the thermal conductivity tensor and a and a are thermo-electric 
tensors. Here a is obtained by the Kubo formula as follows and a is obtained 
via the Onsager relatioE^l as 

1 ^ 

a = -a. (77) 

The expectation value of the heat current 3^ caused by the external electric 
field E is 

J?{k,uj) = 5^a^,(k,a;)E,(k,u;), (78) 

V 

and the thermo-electric tensor a is given by the Kubo formula 

$J,(k,za;,)= rrfre'-^-(T,{jJ(k;r)j:(k = 0)}), (79) 

^0 



We have regarded E and —VT/T as external forces and assumed the Onsager relation 
between the responses to these forces. 
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and 

a^.(k, a;) = - [<^%{k, to + t5) - ^%{k, t5)] . (80) 

Here 

J?ik;T) = e^^j^ik)e-^\ (81) 

and 

f{\i) = jd'xe-'^''f{i^). (82) 

9 Jonson-Mahan Formula: Relation between 
Response Functions 

Let us find the relation between the energy current 

and the charge current 

Introducing a two-time function 

J,(M')=(.)iW^-^*v(0). (85) 
the energy current is written aJ^ 

= - lini ^ - 1;) ^ E ^-(^' ^')- (86) 



f'^i 2 V ^ dt 

Since 



To determine we have chosen one of the equivalent forms of the density £kin 
resuhing in the same value after integration L = J d^ziZkin- Here the equivalent forms of 
the current density is introduced via the integration by parts 

where we have neglected the contribution of the surface. 
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we obtain 

Thus is expressed as 

$^,(k = 0,iwA) = e lim / dTe'''^^F{T,r'), (89) 
^'^^ Jo 

where we have introduced a function 

F(r,r') ^ 5^5^^;,.(ct^(r)cp.(r')j:(k = 0)>, (90) 



P cr 



via the Fourier transform of fl85|) . with t;^; = Px/m. Here we have assumed 
that /3 > r > r' > and dropped Tr from fl69|) . Introducing a function 



is expressed as 



(91) 



(92) 



Employing the Fourier transforno 

= ^ E E ^^n')e*^^"^-'"'^'\ (93) 



20 



Neglecting the vertex correction 



Fir, T')^eJ2Y. M'^U^M^'K'.' (0)cp'.' (0))<, 

pp' aa' 

is factorized as 

P o- 

Since 

Mt')cIM) = -G(p,r'), (4,(T)cp.(0)) = G(p,-r), 

for T, r' > 0, 

Fir, t') ^-eY^Y. ^'^(P' ^')^(P' -^)- 
p 

Via the Fourier transform of the propagator 

^(P'^') = ^E^(P'*^"')^""""'' G(p,-r) = i^G(p,^£„)e'^"^ 

n' n 

F(t, t') in this approximation is shown to have the form of 
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with En and being the fermionic thermal frequency, $^3, is expressed as 

e 



Since 



we obtaiirn 



(94) 



(95) 



(96) 



Thus and $2c relateco by and (ES) for electrons. 

If we regard the Cooper pair as the carriecj of charge and heat, we can 
introduce phenomenologically the current operators as 



and 



where 



j^(x) 



j^(x) 



2m*i . 
1 

2m* 



*^(x)(V^(x)) - (V^"^(x))^(x 



V^^(x))^(x) + ^^(x) fv^(x) 



(97) 

(98) 
(99) 



^t(x) = ^ e-'^^-^Pl vl>(x) = J2 e'^^'^q- 

q q 

These are the GL expressions of the currents [S]. The effective charge and 
mass of Cooper pairs are given as 



e* = 2e, 



m* = 2m. 



(100) 



If we consider the vertex correction, we obtain 

F(r,r') = -5]5]«,G(p,T')G(p,-r)j^, 



where j| is the renormalized current vertex so that (|93p holds in general. 
We have employed 



Jo 



This relation is obtained by Jonson and Mahan: Phys. Rev. B 42, 9350 (1990) and 
also discussed by Mahan: Solid State Physics 51, 81 (1998). 

Since /i = for Cooper pairs, the heat and energy currents are identified, j'^(x) = 
u(x), here. 
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Employing the Cooper-pair currents (I97j) and (1^ instead of the electron 
currents (15^ and fl^ . we obtain 

e* ^ 

TO 

and 

<l>^^(k = 0,iu;A) = (^^'n + + iu;A), (102) 

m 

where 

F(r,rO ^ J]<.(P^(r)P,(r')j-(k = 0)>, (103) 
q 

with being the bosonic thermal frequency and f* = qx/m*. Thus $ 
and are related by fllOip and (11021) for Cooper pairs. 



e 

XX 



10 Ward Identities: Exact Results 




Figure 2: Current vertices for quasi-particles: The left black one is the charge 
current vertex. The right gray one is the heat current vertex. The broken 
line with arrow represents the coupling to the external field. 



First we discuss the derivation [1] of the Ward identity for charge current 
vertex at zero temperature. We introduce the three-point function A^, (/i = 
1,2,3,0), defined by 

kl{x,y,z) = {T{jl{z)i^^{x)i,\{y)}), (104) 

where {A) represents the expectation value of A in the ground state, T is the 
time-ordering operator and z = (z, Zq) with coordinate vector z = {zi,Z2, z^) 
and time Zq. Here ipf{x) and ipliv) ^i-re annihilation and creation operators 
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of t-spin electron. The charge current of electrons obeys the charge- 
conservation law 

E ^Jti^) = 0, (105) 

where 

f,{z) = e^\{z)ij^{z) + e^|(^)V;(^). (106) 
The time ordering of three operators results in the summation of 3! terms as 





- xo)9{xo - 


-yo) 


- {j;{z)ij\{y)ij^{x))9{zo 


- yo)0{yo - 




+ {ij^{x)j;{z)^l{y))e{xo 


- zo)0{zo - 


■yo) 


- {xj;\{y)j;{z)x/;^{x))e{yo 


- zo)9{zo - 




+ {,l;^ix)my)j;iz))eixo 


- yoWivo - 


- Zo) 


- {,l^\iy)i^^ix)j;iz))e{yo 


- xo)9{xo - 


- Zo). 



(107) 



Thus the time-derivative of Af, results in 



_d_ 
dz. 



^Ao{x,y,z) = 6{zo - xo){6{xo - yo){[jo{z),i^ti^)]il^l{y)) 

-9{yo-xo){il^l{y)[foiz),M^)])) 
+6{zo - yo)(^9{xo - yo){M^)[joiz),i^l{y)]) 
-9{yo-xo){[fo{z),i^l{y)W^))) 

+(^{^^'^tiM(y)})- (108) 

Adding the divergence in terms of three coordinate variables and using again 
the time ordering operator T we obtain the four-divergence of A^ as 

E ^^m(^' y^ ^) =<^{ [foi^lM^)] ^\(y)W^o - xo) 



dzu 
=0 



+{T{i;^{x) [f,{z),i;l{y)] })5{zo - yo) 



(109) 
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The last term on the right-hand side vanishes due to the charge-conservation 
law (11051) . Only equal space-time commutation relations are non- vanishing, 

[j-(z),^{(y)]5(zo - 2/0) = ei,\{y)5\z - y), (110) 

and 

[foi.z),ilJ^{x)\5{zo - xq) = -eip^{x)6'^{z - x), (111) 
so that the non-vanishing contribution becomes 

J2 g^K(^^ y,^)=- e{T{i,^ix)i;\iy)})5\z - x) 

11=0 ^ 

+ e{T{i^^{x^{y)])5\z-y). (112) 
Introducing the electron propagator G{x,y) as 

G{x,y) = ~t{T{i;^ix)4iy)}), (113) 
this relation is written into 

3 r. 

^A^(x, y, z) = -ieG{x, y)5\z - x) + ieG{x, y)5\z - y). (114) 

Assuming the translational invariance we set ?/ = and introduce the Fourier 
transform as 

Kl{p,k) = [ d'xe-'P^ [ d'ze-''^{T{j;{z)^^{x)ij\{0)}), (115) 



where the four- momentum is defined as p = (p,Po) and k = (k, fco)- The 
left-hand side of flll4p is evaluated as 

|:AA-(x.0,.) = /|5je-/ ^e-|:,t,.Ayp.*). (116) 

and the right-hand side is transformed as 

J d^xe-'P^ J d'^ze-'^'(^-Gix,0)S\z - x) + G{x,0)5^iz)^ 

= -Gip + k) + Gip), (117) 

where 

G(p,^/.*.e-G(.,0). (118, 
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Therefore we obtain 

3 

J2 KK(P^ k) = eG{p) - eG{p + k). (119) 

IJ.=0 

It should be noted that the factor e represents the charge carried by an 
electron and is automatically taken into account by the commutation relation. 
The vertex function Ff, is introduced as 

A^ip, k) = iG{p) ■ r^(p, k) ■ iG{p + k), (120) 

in accordance with the definiton of the Green function flll3p . Then the Ward 
identity for the charge current vertex is obtained as 

3 

fc) = eG{p)-' - eG{p + k)-\ (121) 

^J.=o 

Since the Fourier transform is introduced as 

px = piXi + P2X2 + p^xz - et, (122) 

where e is the energy and t is the time, xo = t and po = — e (and in the 
same manner ko = —00 with uo being the energy of the external field) in the 
zero-temperatute formalism. 

In the finite-temperature formalism we employ the time-ordering operator 
Tj. and consider the three-point function 

A^(x,i/,z) = {T,{jl{z)^^{x)ij\{y)]), (123) 

where the real time zq and the imaginary time is related by = izo 
and {A) represents the expectation value of A by grand canonical ensemble. 
Taking the charge-conservation law fllOSp into account we obtaiij^ 

E d^^>^ ^) =<^-{ [m,i^,{x)\i^\{y)])5{T^ - r.) 

11=0 ^ 

+ {Tr{i^^{x)[3^,{z),^{y)\])5{T-^-Ty). (124) 

The derivative 9Ag(x, ?/, z)/dTz results in the same form as (|108p where, since d/dr = 
—id/dt, the charge-conservation law becomes 
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Here we have used the relation 

The finite-temperature vertex function is introduced as 

a;{p, k) = [- G{p)] ■ r;{p, k)-[-G{p + k)] , (126) 

in accordance with the definiton of the thermal Green function 

G{x,y) = -{Tr{ij^ix)i,\iy)}). (127) 

Then the resulting Ward identity is the same form as fll2ip in the case of zero 
temperature. For the finite-temperature Ward identity the zeroth component 
of the four-momentum is po = —iSn with fermionic thermal frequency and 
ko = —ioJx with bosonic thermal frequency Ux. 




Figure 3: Current vertices for Cooper pairs. 



The extensioio of this Ward identity to the case of Cooper pairs is 
straightforward. In this note we only consider the local Cooper pair with 
s-wave symmetry in consistent with the local attractive interaction Re- 
placing ^^{x) by ^f(a;) = tpi{x)^f{x) and ^{(2/) by "^^{y) = Vt(l/)^|(2/) in 
(11041) we consider the three-point function as 

M;ix,y,z) = (r{j-(z)vl/(x)vl/t(y)}), (128) 

where \E'(x) and ^'^(i/) are annihilation and creation operators of a Cooper 
pair which has a bosonic character. Using the commutation relation 

[jo^(z), ^\y)]6{zo - yo) = 2e^^{y)6%z - y), (129) 



I have derived the Ward identities for Cooper pairs with s-wave symmetry in 
larXiv:11 08.0815 The Ward identities for anisotropic Cooper pairs have been discussed in 
larXiv:1108."5272l The Ward identities for Cooper pairs have been discussed in comparison 
with those for CDW and SDW in arXiv:1109.1404, 

[Nar] = arXiv:110AQ815, arXiv: 1108.5272, arXi v:1109.1404l 
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and 

[f,{z), ^{x)]5{zo - xo) = -2e^{x)b\z - x), (130) 
the four-divergence of is expressed as 

3 r, 

J2 q-M;{x, y,z) = - 2e{T{^{x)^^{y)})6\z - x) 

+ 2e{T{^ix)^\y)})6\z-y), (131) 

by repeating the same calculations as those for deriving flll2p . Here the 
difference between fermion and boson is handled solely by the time-ordering 
operator T so that the expression of the divergence is common to fermion 
and boson. Introducing the Cooper-pair propagator D{x,y) as 

D{x,y) = -i{T{^{x)^\y)}), (132) 

we obtain the Ward identity 

3 

J2 K^'M^ k) = 2eD{qr' - 2eD{q + k)-\ (133) 

for Cooper pairs where is the counterpart of and D{q) is the Fourier 
transform of i5(x,0) with four-momentum q = (q, go)- It should be noted 
that the factor 2e represents the charge carried by a Cooper pair and is 
automatically taken into account by the commutation relation. 

Although the above derivation for Cooper pairs is formulated at zero 
temperature, fll33p also holds at finite temperature with go being a bosonic 
thermal frequency (go = —i^m)- We are mainly interested in the normal 
metallic phase (T > Tc), the Cooper-pair propagator is a fluctuation propa- 
gator in this case. 

Next we discuss the derivation of the Ward identity for heat current ver- 
tex]^ Here we only consider the local interaction^ll We introduce the three- 
point function defined by 

A«(x,y,^) = {Tr{j^iz)^p^ix)^iy)}). (134) 

The heat current obeys the conservation law 

E|-tfW = 0, (135) 

At=0 ^ 

26 See [Ono] = Ono: Prog. Theor. Phys. 46, 757 (1971). 

■^^ The derivation becomes very simple in the case of the local interaction [Kon]. See 
[Ono] , [Kon] , [Nar] for the derivation in the case of non-local interaction. 
[Kon] = Kontani: Phys. Rev. B 67, 014408 (2003). 



21 



where jo{z) is given by (jSZ]) and is the density of ([T]). The four- divergence 
of at finite temperature becomes 



3 



+ {T,{M^) [j^{z),my)]})6{T, - Ty), (136) 



as f ll24p . Since the interaction among electrons is local, [jo {x) , ipi^{x)] = 
[K,^^{x)] and [j^ (y) , ^l{y)] = [fsT, so that using ^ 

[j^ix),^lj^ix)] = -^^^t(^), [j^iy),4iy)] = -^^/M- (137) 

For simplicity of the representation of the Fourier transform, we discuss the 
zero-temperature case for a while. Employing flllSp the four- divergence at 
zero temperaturj^ is written as 

7^A«(x, z) = — y)6\z - x) + — y)5\z - y), (138) 
^ dz^ dxo dyo 

The Fourier transform of A'^ becomes 

d^ze-'^' J d^xe-'P'"" j d^ye'PyA^ix,y,z) = A'^ip,p-k)i27c)^S\-k-p'+p), 

(139) 

where we have assumed the translational invariance (regarded A^ as a func- 
tion of y — X and z — x) and set A^(x, y, z) = A^(x — y,z — x), and 

A^(j9,p-A;) = j d\x-y)e-'P^''-y^ j d\z-x)e-'^^'-'''^A^^{x-y,z-x). (140) 

At the same time we set G{x,y) = G{x—y) and introduce the Fourier trans- 
form 

^ -G{x -y)= I ^^e^^'^--y\p',G{p'). (141) 



dxo J {2iT 



'^^ Starting from 

A«(x,y,z) = (r{j«(z)V;r(x)4(y)}), 

we obtain 

3 ^ 

J2 ^AQ{x,y,z) ^{T{[j^{z),Mx)]^p\iy)})Sizo-xo) 



+(r{Vt(x)[jo^(z),Vj(2/)]})'5(zo-2/o). 
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Then we obtain 

3 

J2 kXiP^P- = PoG{p) - (Po - h)G{p - k). (142) 

It should be noted that the factor po oi' Po — represents the energy carried 
by an electron and is automatically taken into account by the commutation 
relation. By shifting the four-momentum this relation is converted into the 
Ward identity 

3 

^ A;^r«(p + k,p)= poG{p + k)-' - {po + ko)G{p)-\ (143) 
^J.=o 

for heat current vertex. Here we have used the same relation A^{p + k,p) = 
iG{p + k) ■ ^9(p + k,p) ■ iG{p) as fll20p . This identity also holds at finite 
temperatureo 

The extensionl^ of this Ward identity to the case of Cooper pairs is also 
straightforward. The Ward identity for heat current vertex of Cooper pairs 
is 

3 

+ k,q) = qoD{q + k)-' - (go + ko)D{q)-\ (144) 

fi=0 

where is the counterpart of F^. It should be noted that the factor go or 
go + ko represents the energy carried by a Cooper pair and is automatically 
taken into account by the commutation relation. 

Finally we discuss an application of the Ward identities. If we substitute 
the quasi-particle propagator^ 

G{p)-' = -Po - ep-|, G{p + k)-' = -Po -ko- ^p+i, (145) 

for the full propagators in f ll21|) . we obtain 

3 

J]fc,F^ = e/co + e(ep+.-ep-t), (146) 

/i=0 

As in the case of the Ward identity for charge-current vertex, the relation 

AQ{p + k,p)=[-Gip + k)]-r^{p + k,p)-[-G{p)], 

should be taken into account. 
3" See [Nar]. 

We take the momenta of G{p)^^ and G{p + k)^^ symmetrically. 



23 



so that for = 1, 2, 3 



A" 



(147) 



where we have employed (jlj) by substituting the renormahzed mass of quasi- 
particles m* for m. In the same manner we obtain 



^J.=o 



'1481 



from ffm . Thus for /i = 1, 2, 3 

1 



" 2m 



-Pfi [PO + (PO + fco)] = Vfi {i£n + 



iujx\ 
2 / 



(149) 



These current vertices fll47p and fll49p are consistent with fl94|) and fl96l) of 
the Jonson-Mahan formula. 

If we substitute the Cooper-pair propagator^ 



for the full propagators in (11331) . we obtain 

3 r 

J]A;^A^ = 2eD(0,0)-ie2 (^^ . 



l + e^(q 



k\2 
2/ 



(150) 



SO that for /x = 1, 2, 3 



2) -("1^2 



A' =-4eD(0,0)-'^'j,. 

In the same manner we obtain 



:i5i) 



(152) 



A« = D(0,0)-ieV[go + (go + A:o)] =-2D(0,0)-ieV(^w™ + ^), (153) 

from fll44p for /i = 1, 2, 3. These current vertices (11521) and fll53p are consis- 
tent with fllOip and fll02p of the Jonson-Mahan formula. 



11 Quasi-particle Transport: Relaxation-Time 
Approximation 

Within the relaxation-time approximation we make a microscopic calculation 
equivalent to the Boltzmann transport in the following. 



In general this Ornstein-Zernike form is expected in the long-wavelength limit. We 
take the momenta of D{q)~^ and D{q + k)^^ symmetrically. 
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Figure 4: Quasi-particle transport within relaxation-time approximation: 
The left Feynman diagram describes the charge response to the external 
electric field. The right describes the heat response to the electric field. 



The relaxation-time approximation corresponds to neglecting the vertex 
correction as shown in Fig. 4. In this approximation is given b\j^ 

$:,(k = 0,^a;,) = -^J2J2vlGi^6^)Gi^e^ + ^u,), (154) 

p n 

where the minus sign is the fermion-loop factor. We employ the analytic 
continuation of the thermal propagator G{ien) defined for pure imaginary 
frequencies to the retarded or advanced propagator, G^{e) or G^{e), 

for real frequencies. The analytic continuation results in G^{e) for e„ > 
and (^^(e) for e„ < 0. We have to calculate the discrete summation 

I'ilCOx) = -IY1 Gi^^n)GitEn + tCOx). (156) 

n 

The summation is transformed into the integral^ on the contour C in Fig. [5] 
as 

F{tux)= / —fiz)Giz)Giz + tco,). (157) 
Jc ^7r« 



The footnote for (j93p shows how (|154p results from the Jonson-Mahan formula. Al- 
ternatively 

<i>L(k=o,*c.A) = -i5]5]r^G(*£„)G(z£„ + zc.A)r|, 

p n 

results in (|154D using the Ward identity (|147p . 

34 See, for example, §25 and Problem 7.6 in [FW]. Since tanh(z/2r) = 1 - 2/(z), 



r f/z z 
-r{iujx)= / — tanh— G{z)G{z + iuJx). 
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c 



c 



Co 



c 



C3 



Figure 5: Contours of the integral: These three panels are three different 
characterization of the same complex z-plane. The solid line with an arrow 
represents the contour of the integral. The horizontal and vertical gray lines 
represent the real and imaginary axes of the complex 2;-plane. The dots on 
the imaginary axis represent the fermionic thermal frequencies. The broken 
lines Co and C\ represent the cuts along Imz = and Im^; = —ux. 



The contour C is transformed into Ci, C2, C3 and C4 without encountering 
the singularity!^ of the propagator as shown in Fig. [51 The analytic continu- 
ation of G{ien)G{ien + iuj\) is G^{z)G^{z + iuJx) on Ci, G^{z)G^{z + iujx) on 
C2 and C3 and G^{z)G'^{z + iujx) on C4. Taking into account of the direction 
of four contours we obtain 



Co 



dz 
dz 



f{z)G''iz + toJx) 



G'^(z) - G^(z) 



+ 1 ^f{z)G^iz) G''{z + tux)-G^{z + tux) 



:i58i 



Shifting the variable in the second integral in fll58p we obtain 
de 



I'itoox) 



27Ti 



/(e) 



G''(e + itux) + G''U-iLOx) 



G'^(e)-G^(e) 



(159) 



where we have used the relation /(e — ioox) = /(e) led by e ^'^'^^ = 1 for 
bosonic frequency ux- Since 
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The propagator 



G^(e) -G^(e) =22lmG^(e), 



(160) 



G(z) 



de- 



has a cut along Im 2: = so that G{z)G{z + iuj\) has cuts along Im z = and Im z — —u!\. 
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we finally obtain 

F{uj + i5) = 



°° de 



(161) 



To calculate the DC conductivity we only need the contribution linear in u 

'dG^(e) dG^(e) 



de 

I\u + i5) - r{i5) ■= u I —/(e) 



TT 



de 



de 



ImG^(e) 



2iuj 



de 



d_ 
d~e 



ImG^(e) 



de / dfje] 
vr V de 



ImG^ie 



ImG^(e) 

2 



(162) 



where we have employed the integration by parts. Thus the DC conductivity 
axx is given as 

a,, = hm -[$:,(k = 0,a; + z5) - <l>l,{k = 0,z6)] 



2e' - 



def df{e) 



TT V de 



J2^l[lraG^{ 



1 2 



(163) 



If the life time r is not so short (epT ^ 1 where ep = p'jp/2m), ImG'^(e) 
behaves as —n5{e — ^p) so thalQ 



P/ \ 2 



°° de [■ 1 2 

- ImG^(e) 



(164) 



After the integratioro over e we obtain 



P/ \ 2 



(165) 



See, for example, §8.1.2 in Mahan: Many-Particle Physics, 3rd edition (Kluwer Aca- 
demic/Plenum, New York, 2000). 
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Using the definite integral 



dx 



for a > we obtain 



de 



ImG^^(e) 



_^(a;2 + a2)2 2a3 ' 
l/2r 



(6-ep)2 + (l/2r)2 



TTT. 
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which is equivalent to 
In the same mannei 



as fll54p is given by 



$«(k = 0,^a;A) = -^5^5^t^l(^£n + ^)G(z£„)G(z£„ + ^a;A). (166) 



p n 



The discrete summation 

n 

is transformed into the integral 

""^ de 



(167) 



TT 



-/(e) 



G^(e + u;)+ (e-^)G'^(e-a; 



(168) 



in the same manner as (11611) . The cu-linear contribution is evaluateco as 



de 



de 



de 



1 

+ 2 



A/ a/(€) 



TT V de 



ImG^(e) 



1 2 



G«(e) - G^(e) 
(169) 



Thus 



a^^(k = 0,w-^0) = lim — [$£^(k = 0,a; + i5) -<l>i.(k = 0,z5)] 



2e 



rfe/_9/(e) 



ImG^fe 



(170) 
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Using (|TTfl) and ([Tig| 



(k = 0, = E E r?G(ze„)G(*£„ + iLoxWl, 



results in (I166p . 

Using (|160p . the right-hand side of the first hne of (I169P is equal to 



de 



/(e) 



2e4-ImG^(e)+ImG^(e) 



ImG' 



{e)=^u: j" [|(/(e)..)] +/(.)} [lmG«(6) 



where we have employed the integration by parts. 
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In the same manner as (11641) we obtai: 



a. 



(k = 0,a;^0) ~ 2e^ 



(171) 



which is equivalent to fj^T]) . 

As another way of evaluating the summation in fll56p we employ the 
spectral representation fl20|) and obtair0 

p{y)dy 



isn + iuj\-y 



dx 



dx 

—oo J —oo 



dy- 



^^^{x)p{y^ly 
X — y + iu\ B ^-^ 

n 

p{x)p{y) 



1 



iSn-x iSn + iuJx-y 



[f{y)-f{x)\. 



(172) 



x-y + iujx 
The imaginary part is 

/oo /"OO 
dx I dyp{x)p{y)[f{y) - f{x)]6{x - y + u) 
-oo J — oo 
oo 

dxp{x)p{x + u) [f{x + u) — f{x)~\ 
dxp{x)[p{x — oj) — p{x + u)~\f{x). (173) 
The w-linear contribution is estimated as 



— TT 



-TT 



'OO 

oo 



/oo 
dxp{ 
-oo 



dp{x) 



:i74) 



Employing the integration by parts we obtain the final result equivalent to 
(HSil) 

df{x) 



im[r{u + iS)] ■= 



TXUO 



dx 



p{x) 



dx 



(175) 



40 Pornially axx is proportional to the charge of electrons. However, the sign of the 
summation at low temperature is determined by that of A^i where the density of states 
N{^p) is expanded as N{£_p) '—. N{0) + Ni ■ ^p. In the presence of particle-hole symmetry, 
A^i and thus axx vanish. See (10.19) in 2 for the discussion on the symmetry. 

''^ Here the relation 



S-^Q B ^ I. 



is employed taking into account that e„ + wa is also a fermionic thermal frequency. See, 
for example, §25 in [FW]. 
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12 GL Transport: Gaussian Fluctuation 

Here we discuss the GL transport jfl We consider the GL free-energy density 



2m* 



V ' 



(176) 



in the Gaussian approximation. Introducing the Fourier transform of the 
complex order-parameter field 

v^(x) = J2 e^'^'^q' ^* W = E e""'"^;, (177) 

q q 

the GL free-energy E is given by 

E = J d?xE{^) = i: (« + ^) (178) 
q 

The Gaussian fluctuation is evaluated a^ 

T T 1 
/vl/*vl/ \ = = flTQ) 

\ ^ a + q'/2m* a 1 + ^ V 

where (A) is the expectation value of A under the Boltzmann weight exp(— E'/T). 
Here the correlation length ^ of the fluctuation is introduced as 

= irV- (180) 

The relaxation of the fluctuation is determined by the time-dependent GL 
equation 



q 



In terms of the relaxation time Tq 



- = -(« + /7)=-(l + eV), (182) 
Tq 7 V 2m* 77^ ' 



^"^ We follow the discussion for the conductivity by Skocpol and Tinkham: Rep. Prog. 
Phys. 38, 1049 (1975) and Tinkham: Introduction to Superconductivity, 2nd edition 
(McGraw-Hill, New York, 1996). Then we extend the discussion to the case of the thermo- 
electric tensor and our result (|197l) is identical to that in [3] and [S]. 
The Gauss integral for complex variable z is performed as 

dz*dz _ f dwfo („2^^2) _ 1 

27n ~ ' 
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(11811) is written as 



1 



dt 



The Kubo formula for classical variables is 

1 



T 



j:{o)j:{t))dt, 



:i83) 



:i84) 



where the charge currentcj is given by 



:i85) 



The current- current correlation 



is evaluated by using 



:i86) 



-2t/rq 



:i87) 



where the time-dependence of the order-parameter field is determined by 
as 

^q(t) = ^qe-*/^'^. (188) 



After the time-integral we obtain 

e* \ 2 1 



;i89) 



44 



and 



The charge and heat currents corresponding to (jHT]) and (j98l) are given by 



1 



2m* 



^(V**(x))4'(x) + ^'*(x)(v^'(x)' 
as (4) and (5) in Using (|177l) the Fourier-transformed charge current is 



m* ^-^ V 



as (j73p . Thus the x-component of the uniform current J| = ]% (k = 0) is given by (jl85p . 
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Using ( I179p we obtain the DC conductivity 



In 2D the q-summation is performed ai 



d^^^^^ = T7^' (191) 



q 

so that 



T 

"""" " 16^T^- ^^^^^ 



- c 

Next we consider the Kubo formula for the thermo-electric tensor a 



where the heat current is given by 

d d 



1 1"°° 

- j (J'MJS(t))dt, (193) 



q 

which corresponds to (jHSD- Although in fll88p Tq is assumed to be real 



here it is extended to be complejo and using its complex-conjugate r* the 
Here we consider the film whose thickness is d and assume d ^ ^ so that 

The replacement — !■ q'^ /2 can be done for the present symmetric case in 2D. Since 



we obtain 



{x + 1)3 2 



^® If Tq is real, axx vanishes. Since holes carry opposite charge to electrons, axx and axy 
that are odd in the electric charge vanish by the cancelation between contributions from 
particle- and hole- charge currents in the presence of the particle-hole symmetry. On the 
other hand, such a cancelation does not work for the heat current. Even in the presence of 
the particle-hole symmetry axx and axy that are even in the electric charge do not vanish. 
See, for example, Niven and Smith: Phys. Rev. B 66, 214505 (2002) for the consideration 
in terms of the symmetry. 

In general the complex GL equation 

^ = * + (1 ici)V2§ -{1 + iC2)\'i'\^^, 
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time-dependence of the conjugate of fll88p is given as 

^*(t) = ^*e-*/"^. (195) 

Since 

we obtaii]|f§ within the hnear ordeiQ of 7"/7' 



where we have assumed 7' ^ \^'\ with 7 = 7' + ^7" and used f llSOp . 



J50 



In 2D the q-summation is performed a; 
so thatEll 

kl 7" 7a 

ttxx — . r In — ;^r, (199) 

where we have used the Onsager relation fl77|) and introduced Ta as (^cO^ = 
Ta/(T-T,). 
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Figure 6: AL process of Cooper-pair transport. 



13 Cooper-Pair Transport: AL Process 

Here we derive the GL results (11921) and (11991) from the AL process of the 
microscopic calculation. Such a process corresponds to the relaxation-time 
approximation neglecting the vertex correction. Since the effective interac- 



is employed to study the dynamics of the complex order parameter 5*. The origins of the 
imaginary part of Tq are discussed in §9.3 of [3j for example. The non-vanishing imaginary 
part is related to the violation of the particle-hole symmetry. 

The Drude-like formula (|189l) was extended to the case of the thermo-electric tensor 
by Howson, Salamon, Friedmann, Rice and Ginsberg: Phys. Rev. B 41, 300 (1990) and 
Ausloos, Clippe and Patapis: Phys. Rev. B 46, 5763 (1992). However, their discussion is 
incorrect so that they could not reach (jl97p . 

'^^ We have put = except the numerator of the right-hand side of (|196l) . Thus the 
time-integral is identical to that leads to (I189p . 
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tion L also obeys the same 





51 Our resuh is identical to (4.31) in [3] and TABLE I in [S]. 

5^ The Ward identities for L are given as 



3 



k^Kliq, k) = 2eL{q)-' - 2eL{q + k)-\ 



and 



3 



5^fc^A^(g + fc,(j) 



qoL{q + k) 



1 



(go + ko)L{q) 



-1 



and we employ 




as (fT50)) . 
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current verticeq^ for L as 

= AeNmoQ,. (200) 

and 

A« = 2NiO)eoq,[toJrr. + (201) 

These relations also hold even when we take the effect of impurity scattering 
into account 1^ 

The charge response to electric field is determined by 

<i>l,{k = 0,iux) = ^5^5^A^L(^a;™)L(^a;^ + ^a;,)A:, (202) 

q m 

where we have used the simplified notation L{iuJm) = L{q^, ium)- The discrete 
summation ^ 

P{iu}x) = ^ L{iujm)L{iujm + (203) 

m 

is transformed into the integral^ 

P{iujx)= / —n{z)L{z)L{z + iuJx), (204) 
Jc 

where the contour C is shown in Fig. 5 and 

n{z) = -J-^, (205) 

is the Bose distribution function. Employing 

L^{x) - L^{x) = 2i ImL^ix), (206) 
and the relation n{z — iuj\) = n{z) led by e"*^"^^ = 1 we obtain 

/°° rlr r i 

—n{x) L^(x + u) + L^{x - u) ImL^(a;), (207) 
-oo L J 



The frequency factor in (|20ip appeared in (34) of [Uss]. However, such a factor 
cannot be derived from the preceding discussion in [Uss]. I shall discuss this point in the 
supplement noticed in the footnote 63. 

[Uss] EE Ussishkin: Phys. Rev. B 68, 024517 (2003). 

See, for example, §8.3 in j^. The value of becomes a function of the mean free 
path due to impurity scattering. 

55 See, for example, §25 in [FW]. Since coth(z/2T) = 1 + 2n(z), 



r fiz z 
r{iujx)^ —coih—L{z)L{z + iuJx)- 
J c 
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by repeating the procedure leading to (11611) . The contribution hnear in u is 
obtained as 



where 



dx ( dn{x) 

TT 



1 



dx 

TqX 



ImL^ix) 



A^(0) + {r^xf 



(208) 
(209) 



with ?7 = e + ■ Thus the resuhing DC conductivity is 



2 /"°° dx ( dn{x) 



A 



dx 



ImL^ix 



AttT 



/■oo 

/ dx 


sinh — - 


-2 


ImL^(x) 


1 — oo 


L 2tJ 







(210) 



In order to discuss the effect of the low-energy critical fluctuation at finite 
temperature it is enough to use the high-temperature expansioirflsinh(x/2T) = 
x/2T under the assumption e <^ 1. By performing the integraPl we obtain 



c^xx = SeVoT ^ 



4^2 

X 



In 2D the q-summation is performed ar^®' 



E 



e + eo'q2)3 Svrrf J, {x + e) 



dx- 



X 



1 1 



IGnd e 



(211) 



(212) 



The high-temperature expansion is compatible with the Kubo formula (I184|) and 
P93p for classical variables. 



" Since 



sinh — - 
2T 



the integral is evaluated as 



■"^"(^'l ■=■ (^) ( 



da; 



sinh ■ 



2T 



ImL^(x) 



iV(0)/ V 772 + (Toa;)2 
2T \2 

2to 



7V(0) 



(y^ + ?7^) 



2^2 ■ 



Using the definite integral appeared in the footnote for (|165p we obtain 

dy TT 1 



By the same procedure as the footnote for p9ip 



dx 



1 



(a; + e)3 2e' 
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so tha' 



59 



e2 roT e2 1 



27id e 16d e 
The heat response to electric field is determined by 



(213) 



$i(k = 0,^a;A) = ^5^5^A«L(za;^)L(^a;^ + ^a;A)A:. (214) 



q m 



The discrete summation 



I^{iuJx) = ^ XI + ^-^^L{iu^)L{iUm + «wa), (215) 



is transformed into the integral 



-nix) 



n 



x + '-^]L^(x + uj) + {x-'-^]L''(x-u 



In the same manner as (11691) we obtain 

I^iu + ^6)-I^i^6)=^u H^f^^ 

. TT V dx 



X 



ImL^(x) 



ImL^(x). 
(216) 

(217) 



Employing the high-temperature expansion n{x) '= T/x, (12171) reduces to 



47r 



dx 



L^(x) - L'^(x) 



IX 



Here 



mx) 



L^(x) 



(218) 



(219) 



N{0)r]-iTox' ' ' N{0)r] + iT^x' 

with To = Ti+iT2 and Tq = Ti —iT2 and we can put Ti ^ \t2\. The right-hand 



Here we have used tq = vr/ST derived within the ladder approximation. See, for 
example, §6.2 in [3]. 
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side of ( 12181) is evaluated exactly by the residue sJ^^' 

T T2 1 



— lUJ- 



2iV(0)2rir/2- 

Thus the analytic continuation of f l214p is given af^^^' 



(220) 



C(k = 0,a; + z5) = Ha;)4eT^5^ 



In 2D the q-summation is performed as (11981) : 



(221) 



2„2\2 



8nd 



dx- 



X 

(x + e) 



1 , Xc 
: m — , 

87rci e ' 



so that 



2nd Ti e 
By the Onsager relation (1771) we finally obtain 



_ 1 ^ 



m 



2Tid Ti T — Tc 



(222) 



(223) 



(224) 



which is equivalent to (I199p where T2/T1 = ^' and Xc is chosen as Xc/e = 
Ta/(T-T,). 



60 



The integral 



dx 



L^(x)-L^(x) 



L^{x)~L^{x) 



IX 



is written as 



^0 +^"0 



dx 



1 



1 



iV(0)2 

and is evaluated by the residue as 



rj — itqx rj + iTqX 



(?7 - iTox){ri + iT*x) ' 



2tT Tq* - Tq 1 

iV(0)2r*+To 772 ' 



where the position of the pole is determined by the condition ti 3> |t2|. 

Since T2 < in 3D for our free electron dispersion (U), (|22ip implies that a^x is 
proportional to the charge of Cooper pairs. For example, (2.27) in Fukuyama, Ebisawa 
and Tsuzuki: Prog. Thoer. Phys. 46, 1028 (1971) shows T2 < after the correction 
pointed out by several authors: its right-hand side should be multiplied by —1. 
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14 Ward Identities: Ladder Approximation 



It is instructive to show how the Ward identity is satisfied in the ladder 
approximation that is usually employed to obtain the explicit form of the 
Cooper-pair propagator. In particular it is very important to clarify how the 
Ward identity for heat current vertex is satisfied, because it will clear away 
some confusion^ seen in the literatures. Such a clarification becomes in- 
evitably complicated so that it will be uploaded as a separate supplementary 



15 Current Vertices: Perturbational Analy- 
sis 

As has been shown in this note we do not have to calculate the current 
vertices for Cooper pairs to obtain a^x and a^x- However, it is instructive 



included in the supplementary note noticed in the previous section. 

16 Exercise 

Find and correct the errors in this note. (Do not trust the results in this 
note before you check them.) 

For example, if we employ the heat current vertex for electrons with the frequency 
factor i£„+iWA/2 in perturbation expansion, we obtain some wrong results. The frequency 
factor should be accompanied by the full propagator. It should not be accompanied by 
the free propagator. 

^■^ I shall upload a note entitled A Diagrammer's Note on Superconducting Fluctuation 
Transport for Beginners: Supplement besides the three notes in series. 

If we perform the perturbational calculation of the current vertex without the ground 
as the Ward identity, it is difficult to convince ourselves of the validness of the result. For 
example, there are many errors in the calculation of vertices by Varlamov and Livanov: 
Sov. Phys. JETP 71, 325 (1990) and 72, 1016 (1991), besides the fatal error, concerning 
the "cancelation" , corrected by [RS] . 

[RS] = Reizer and Sergeev: Phys. Rev. B 50, 9344 (1994). 

The factor of the heat current vertex (45) in [RS] is incorrect but this result is 
supported by Serbyn, Skvortsov, Varlamov and Galitski: Phys. Rev. Lett. 102, 067001 
(2009). I do not understand the reason of the support. On the other hand, the heat 
current vertex (21) in [Uss], which is cited as (10.23) in [5] and corresponds to (|20ip with 
ujx = 0, is insufficient to obtain axx- 



note. 
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